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A systematic construction of the Green's matrix for a second order, self-adjoint matrix 
differential operator from the linearly independent solutions of the corresponding homoge- 
neous differential equation set is carried out. We follow the general approach of extracting 
the Green's matrix from the Green's matrix of the corresponding first order system. This 
construction is required in the cases where the differential equation set cannot be turned to 
an algebraic equation set via transform techniques. 

I. INTRODUCTION 

In physics, matrix differential operators acting on vector functions appear in many different con- 
texts from classical electromagnetism to quantum field theory. Green's matrices of these operators 
are needed because of their own physical interpretation as propagators in quantum field theory, or 
in order to find the solutions of the corresponding non-homogeneous differential equation set. 

In most of the cases. Green's matrices are obtained by using the Fourier transform technique 
or by using eigenfunction expansions which turn the differential equation set to an algebraic one. 
However, these techniques are not applicable in some circumstances such as for the matrix dif- 
ferential operator appearing in the 1+1 dimensional Abelian-Higgs model [1, 2]. In this model, 
when small field fluctuations around classical field configurations are investigated, Lagrangian of 
the theory, which is second order in fluctuations, involve a 4 x 4 matrix differential operator. Di- 
agonal entries of this operator are differential operators of the modified Bessel type. Off-diagonal 
"potential" terms are functions of classical field configurations which are only available as discrete 
numeric data for generic values of parameters in the theory. Green's matrix of this operator is 
required in calculating the functional determinant of the operator which gives one-loop corrections 
about a classical solution such as an instanton (For a nice account of functional determinants using 
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Gel'fand-Yaglom technique see [3]). Green's function technique is one of the methods used in such 
a determinant calculation. Due to the existence of the modified Bessel type operators and the 
discrete numerical data, it is not possible to apply Fourier transform and eigenfunction expansion 
in this case. However, in a numerical study, it is relatively easy to obtain linearly independent 
solutions of the corresponding homogeneous differential equation. Therefore, construction of the 
Green's matrix from these solutions as in the case of the single differential operator is required. 

It is worth considering the underlying physical problem in some detail in order to see why one 
encounters a matrix differential operator which is hard to handle. In [4], 't Hooft studied the one- 
loop tunneling amplitude in the background of a Yang-Mills instanton for a theory which contains 
massless scalar and fermion fields . In this calculation, field fluctuations do not couple, and the 
functional determinants of single differential operators are calculated. Dunne et al. [5, 6] extended 
the instanton determinant calculation to the arbitrary quark mass case, 't Hooft [4] pointed out 
that in order to remove the infrared divergence of the theory, one needs to introduce the Higgs 
field. However, due to simple scaling arguments, there is no instanton solution in this case. It is 
still viable to do the calculations in the same instanton background, but with certain instanton size, 
since it was shown in [4] and in [7], with a more elaborate discussion, that the Higgs particle can 
be taken as approximately massless. However, if one studies the effect of the quantum fluctuations 
around an instanton background with non-trivial Higgs field configuration, as in the case of the 
1-1-1 Abelian-Higgs model, field fluctuations do couple to each other and one needs to struggle with 
the functional determinant calculation of a matrix differential operator. 

In the physics literature, a general construct/ion of the Green's matrix from the set of solutions of 
the corresponding homogeneous differential equations does not seem to exist. Among the standard 
references of mathematical physics, only Courant and Hilbert [8] involves a short discussion on the 
properties of the Green's matrix ("tensor" as called in [8]) without a construction. Baacke [9] gave 
a construction for a specific matrix differential operator in a heuristic way. He studied one-loop 
effects in various field theories using Green's matrices found by this construction [1, 2, 9, 10]. In 
[11], again there is a construction for a specific operator with the main emphasis on the boundary 
conditions of the underlying physical system which is a magnetic multilayer structure. 

In the mathematics literature, construction of Green's matrix for a second order, self-adjoint 
differential operator from the solutions of the corresponding homogeneous differential equation 
set does exist. Naimark [12] studied Green's matrices of the -n}^ order, linear matrix differential 
operators for general homogeneous boundary conditions relating vector function and its derivatives 
up to the (n — l)**^ order at the boundaries. He gave the Green's matrix form for this general 
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system and outlines a way to prove his result. These results and analysis are too general, so a 
special study of the physically relevant case, which is the second order and self-adjoint operator, 
is still valuable. Let us mention several other related works. Bhagat [13, 14] worked on the case 
of second order, self-adjoint, 2x2 matrix differential operator. Heimes [15] gave an analog of our 
result for second order, linear systems without any expUcit construction, and only mentioned that 
the proper method is to transform the second order system to first order system. Jodar [16] worked 
on an algebraic construction which may not work on every case. 

In this paper, we consider a generic second order, self-adjoint [20] n x n matrix differential 
operator of the form 



Mil,,. Vuix) ... Vmix) 

Vuix) M22,x ••• V2nix) 



Vl„{x) V2n{x) 

where the diagonal entries are of the form 



nn,x 



(1) 



y- (pi (a;) y- ) + Qi (x) 
ax \ ax 



and the off-diagonal entries Vij = Vji are just continuous functions. First of all, we develop the 
properties of the Green's matrix in Section II. Section III is devoted to the construction of the 
Green's matrix. Construction is carried out in two ways. In the first way. Green's matrix of second 
order system is extracted from the Green's matrix of the corresponding first order system. In this 
construction, general approach developed in Cole [17] is followed [21]. This construction relies on 
the same basic idea as the constructions of [11] and [15]. In the second way of construction, we 
start with a guess on the form of the Green's matrix, which is along the line of [9]. 

To fix the notation, let us note that small bold letters, e.g. y, represent vectors; capital bold 
letters, e.g. G represent matrices, x appearing as an index refers to a differential operator such 
as M-r. Repeated indices on different matrices are to be summed over, unless otherwise stated. 
Repeated indices on a single matrix refer to a diagonal element. And a word about the nomenclature: 
to distinguish the Green's function of a single differential operator, we choose the name Green's 
matrix for coupled equations. 
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II. PROPERTIES OF THE GREEN'S MATRIX 

Let's consider a linear, second order, coupled differential equation set; 

M^y{x) = h{x), xe[a,b], (2) 

where y , h are n dimensional vector functions, and M^; is n x n dimensional, self-adjoint matrix 
differential operator of the form (1). Green's matrix, G{x,t) of the differential operator, M^; can 
be defined with the formal solution 



y(x) = / dtGix,t)h{t), 

J a 



(3) 



where G {x,t) is n x n matrix. 

In this paper, homogeneous boundary conditions are considered: 

y (a) =o, y (6) = o, 

where o is the n dimensional zero vector. These boundary conditions impose the following conditions 
on the Green's matrix; 

G(a,i) = 0, G{b,t) = 0, 

where O is ra x n zero matrix. 

The formal solution of the differential equation set implies 

Mx,ijGjk {x, t) = SikS {x-t), (4) 

where i, j, k indices run from 1 to n. This relation is an equality of distributions. Integrals of these 
distributions with a test function yield further properties of the Green's matrix. In obtaining these 
properties, i = k and i k cases will be investigated separately. 
For i ^ k: 

k^^ column of the Green's matrix is a solution of the homogeneous differential equations except 
the k^^ equation, as implied by 

Mx,ijGjk {x, t) = 0. 



Each differential equation involves a term of the form 



Mx,iiGik {x,t) 



Gik {x,t) 
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where there is no summation on i. In order to satisfy these (n — 1) homogeneous differential 
equations, first and second order derivatives in the above term should not yield any singularities, 
since the other terms of the differential equation contain continuous potentials. Thus, elements of 
the k^^ column should be continuous and should have continuous first derivatives for any x € [a, b], 
except Gkk {x, t) which is investigated m i = k case below. 
For i = k: 

k^^ column of the Green's matrix is a solution of the k^^ homogeneous differential equation for 
X G [a, t) U {t, b] , as implied by 

Mx,kjGjk (x, t) = d{x-t), 
where there is no summation on k. This equation contains the term 

Mx,kkGkk {x, t)= (^pk (x) + Qk (x) Gkk {x, t) . 

Since 

Mx,kjGjk {x, t) = 0, 

for X ^t, Gkk ix,t) should be continuous and has continuous first derivatives for points other than 
x = t. 

Let's consider the behavior at x = t. Since all elements of the k^^ column other than Gkk {x,t) 
are continuous at x = t, Dirac delta behavior comes from Gkk {x,t). A discontinuity in Gkk {x,t) 
yields a more severe singularity than Dirac delta upon taking the second derivative. Thus, Gkk {x, t) 
is continuous at x = t and the first derivative of Gkk {x, t) has the usual discontinuity 



lim -^Gii {x,t) 
e->o dx 



t+e 



t-e Pi (*) ' 

As in the case of the Green's function for a single differential operator, self-adjointness of Mj; 
and the homogeneous boundary conditions yield a symmetry property for the Green's matrix: 

(^'^2) Mx,ijGjk {x,xi) = (G^);. {x,X2) 5ikd {x - xi) , 

^ Gil (x, X2) Mx,ijGjk {x, xi) = Gki {x, X2) 5{x- xi) , 

and 

^ Gik {x, xi) Mx,ijGji {x, X2) = Gik (x, xi) S{x- X2) ■ 



After integrating the above two equations over the interval [a, b] and subtracting them side by side, 
one obtains 

/ [Gil {x,X2) Mx,ijGjk {x,xi) - Oik {x,xi) M^^ijGji {x,X2)] dx = Gm {xi,X2) - Gik {x2,xi) . 

J a 

Calculating the left-hand side: 
For i = j: 

fb 



Gil {x, X2) (pi (x) -^Gik {x, xiU - Gik {x, xi) (pi (x) -^Gii (x, X2] 



dx. 



Adding and subtracting Pi (x) ^Gu {x, X2) 4-Gik {x, xi) yield 



^ (gu {x, X2)pi (x) ^Gik (x, xi) j - ^ \ Gik (x, xi)pi (x) ^Gu (x, X2] 



dx. 



After the integration, one obtains 



Gii{x,X2)pi{x) ^Gik{x,xi) \ - \ Gik{x,xi)pi{x)-^Gii{x,X2) 



x=b 



0, 



from Gjk (a, x') = 0, Gj^ (6, x') = for any j and k. 
For i ^ ?": 

rb 

X] / i^ii ^2) Mx,ijGjk {x, xi) - Gik {x, xi) M^^ijGji (x, X2)] dx, 

contains terms like 

I J [Gil {x, X2) Mx,i2G2k {x, xi) - Gik {x, xi) Mx,i2G2i {x, X2)] dxj + 

I j [G21 {x, X2) Mx,2iGik (x, xi) - G2k {x, xi) Mx,2iGii (x, X2)] cixj , 

which vanish, since the matrix differential operator, M-^ is symmetric. 
Thus, one obtains the symmetry property of the Green's matrix: 

Gkl {xi,X2) = Gik {x2,xi) . 

As a result. Green's matrix of a second order self-adjoint matrix differential operator satisfies 
(4) and the homogeneous boundary conditions. Rewriting them together, we have: 

Mx,ijGjk (x, x') = 6ikS (x - x') ; Gjk (a, x') = 0, Gjk {b, x') = 0. 



Properties of the Green's matrix developed in this section can be summarized as: 
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k^^ column of the Green's matrix satisfies the homogeneous differential equations except at 
one point x = t for equation i = k. 



Mx,ijGjk (x, t) = 0, X £ 



[a, t) U{t,b], i = k, 
[a, b] , i ^ k. 



• Green's matrix is continuous at x = t. 

• Derivative of the Green's matrix at point x = t is continuous for the off-diagonal elements 
and has a jump of 1/pi (t) for diagonal elements. 



lim -^Gii (x, t) 



t+e 



t-e 



i^j 



• Green's matrix has the the following symmetry: 



G{x,t) = {t,x). 



These properties are also given in [8]. 



III. CONSTRUCTION OF GREEN'S MATRIX 



A standard way of constructing the Green's function for a second order, linear, self-adjoint 
differential operator 

^^"^(^^")^)+^^")' 

defined in [a, b] is to use the two linearly independent solutions of homogeneous differential equation 
satisfying, 

LxU (x) =0, u (a) = 0, 
L^v (x) =0, V (b) = 0. 

The motivation for such a construction follows the observation of two points. First, Green's function 
satisfies the homogeneous differential equation, except at x = t. Second point is the correspondence 
between the derivative property of Green's function; 
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and the Wronskian of the solutions u and v; 

/ / constant 
W [u, V) = uv — vu = . (5) 

Since the columns of our Green's matrix satisfy the homogeneous differential equation set except 
at one point x = t, it is suggestive that the Green's matrix can be constructed from the solutions 
of the homogeneous differential equation. In this section, this construction will be given. In 
Section III A, analogs of the Wronskian, (5), are obtained. A direct approach for constructing the 
Green's matrix involves first transforming the second order differential equation set to a first order 
differential equation set. Then, Green's matrix of the second order set is extracted from the Green's 
matrix of the first order set. This approach is handled in Section IIIB. 

A. Analogs of the Wronskian 

In general, 2n linearly independent solutions of 

M^y (x) = o, 

can be (re)defined in such a way that n of them satisfy the left boundary condition, and the others 
satisfy the right boundary condition. Let's call them as u" and v^, respectively, satisfying 

u" (a) = o, v'' (6) = o, 

where Greek superscripts label the solutions. 

Following the similar steps leading to (5), it is possible to obtain analog relations for the matrix 
differential operator case. One can write 

< (x) M,,,jv^ (x) = 0, 
v^{x)M,,iju'^{x) = 0. 

Subtracting side by side yields (a and /3 superscripts are suppressed since the equation holds for 
every a and p. Also, x dependence of the solutions are not explicitly shown up until the final 
result.) 



UiMx,ijVj - ViMx,ijUj = 0. 
Since the matrix differential operator is symmetric, terms like U1M12V2 and V2M21U1 cancel each 
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other. After these cancellations, one obtains 

^ {uiMx^iiVi - ViM^^iiUi) = O^^^Ui {piv'i}' - Vi [piu'i)' = 0, 

i i 

i 

=> (uiPiv'i — ViPiu'j) = constant. 

i 

After putting the superscripts which label the solutions, and showing the explicit x dependence 
of solutions, one ends up with 

Y,Vi {x) (u- {x) j-J {x) - {x) ^uf (x)] = C-^ (6) 

where C"^ are constants. Matrix form of this equation is 

(x) P (x) V' {x) - {\J'f (x) P (x) V (x) = C, (7) 

where U and V are n x n matrices whose columns are u"^ and vectors, respectively; and P 
matrix is defined as 

P (x) = diag [pi (x) ,p2 (x),... ,pn (x)] . 

Other two relations that can be derived similarly are, 

(x) P (x) U' (x) - (U')^ (x) P (x) U (x) = O, (8) 
(x) P (x) V (x) - (V)'' (x) P (x) V (x) = O. (9) 

Here, using the boundary conditions u" (a) = o and v'^ (6) = o in (8) and (9), respectively, one 
finds that the constant matrices on the right-hand sides equal to zero. Rearranging these equations 
yields symmetric matrices 

PU'U"^ = (U'U-^P)^ , (10) 
PV'V-^ = (V'V-ip)'', (11) 

in (a, 6), since the P matrix is diagonal. Using these symmetric forms in (7) yields 

PV'V"^ - PU'U"^ = (U^) Cy-\ (12) 

in (a, 6). Note that right-hand side is necessarily a symmetric matrix due to symmetry of the 
left-hand side. 
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B. Green's matrix of the second order operator from the Green's matrix of the 
corresponding first order differential equation set 



Let's rewrite the second order system given in (2) more explicitly 

Mii^^yi [x) + 1^12 (x) y2{x) + ... + Vin (x) y„ (x) = hi (x) , 

Vl2 (x) yi (x) + M22,xy2 (x) + . . . + V2n (x) 2/„ (x) = /l2 (x) , 

: + : + : + : = : 

Vln (x) yi (x) + V2n (x) ^2 (x) + . . . + Mnv^xVn {x) = K (x) . 
with the boundary conditions 

yi (a) = 0, yi (b) = 0. 

Let's rewrite the system of second order differential equations as a system of first order equations 
with the definitions 

Zi (x) = yi (x) , Zn+i (x) = ^Vi (x) . 
Then, the first order system can be put in the following form; 

z' (x) = A (x) z (x) + f (x) 

where 



A = 






.. 


. 


1 


.. 


. 




Zl 










.. 


. 





1 .. 


. 




Z2 









— 31 

PL 


.. 

_Vl2 

pi 




Vln 
PI 




-A 

pi 


.. 
.. 


. 1 

. 


, z = 


Zfi 
Zn+1 


,f = 




hi 

Pi 


_V21 
P2 




Vln 
Pi 





Pi ' ' 


. 




Zn+2 




hi 

P2 


Vni 

Pn 


Vn2 
Pn 


In 
Pn 





.. 


Pn . 




Z2n 




h„ 

Pn 



The boundary conditions can be restated in this notation 

BaZ (a) + BfoZ (b) = o, 
where Bq and Bf, are 2n x 2n matrices in the following form 





I 


O 




O 


o 


B„ = 






, Bfo — 








O 


O 




I 


o 



(13) 



(14) 



(15) 



where I is the n x n identity matrix. 

Green's matrix of the first order system can be defined with the formal solution 



= f dtGi{x,t)f{t). 

J a 



Then, the formal solution of the second order system is 

rb 



y (x) = U„x2nZ (a;) = / dt'Unx2nGi{x,t)f (t) 
J a 



where 



nx2n 



I o 



f (t) can be rewritten as 



f (t) = L2nxnP ^ {t) h {t) , L2nxn = 



O 
I 



Then, 



y{x)= / dt\Jnx2nGi{x,t) 
J a 

Comparing this result with (3) yields 

G {X, t) = Unx2nGi (X, t) L2nxnP~^ {t) . 
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(16) 



Multiplications with U„x2n and ^2nxn choose up-left nxn block of the Green's matrix of the first 
order system. 

The relation between the Green's matrix of the first order system and the Green's matrix of the 
second order system is established. Let's continue with reproducing the result of Cole [17] for the 
Green's matrix of a first order system. 



1. Green's matrix of a first order differential equation set 

Let's have a "generic" first order differential equation set in the form of; 

z' (x) = A (x) z {x) + f (x) , 

where z, f are / dimensional vectors and A is an / x / dimensional matrix. 
First, assume a particular solution of the form 

Zp (x) = W (x) g (x) 



(17) 



12 

where g is an unknown column vector, and the so called fundamental matrix, W is the matrix 
whose columns are the / linearly independent solutions of the homogeneous differential equation 
set, z' (x) = A (x) z (x); i.e. 

W (x) = A (x) W {x) . 

Then, putting the guess for the particular solution in the non-homogeneous equation yields the 
formal solution for g (x) as 

g(x)= r dtw-Ht)Ht)- 

J a 

Therefore, the general solution for the generic first order differential equation set is 

z{x) = W (x) / dt W"^ (t) f (0 + W (x) c, (18) 

J a 

where c is a constant vector. 

Any well-posed boundary condition which yields unique solution to the boundary value problem 
can be put in a matrix form. Consider the boundary conditions 

Baz (a) + BbZ (5) = o. (19) 

Applying these boundary conditions to the general solution fixes c as 

c = -D-^BbW(6) / dtW-^{t)f{t) 

J a 

where D is defined by 

D = B„W(a) + B6W(6). 

Using this result in the general solution, one obtained 

/•X pb 

z (x) = W (x) / dt W"^ it) f (t) - W (x) D'^BfeW (6) / dt W"^ {t) f (t) . (20) 

J a J a 

This result can be put in a final form by rearranging the first term on the right as; 

W (x) r dt (t) f (t) = W (x) Ti-^T) j dtW-^ (t) f (t) , 

J a J a 

= W (x) {B„W (a) + BfcW (5)} / dt W"i (t) f {t) , 

J a 

= W (x) D-^B„W (a) r dt W"^ (t) f (t) 

J a 

nx 

+ W (x) D-^BfeW (6) / dtW-^(i)f (i). 

J a 
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Using this in (20) yields 

z (a;) = W (x) D-^B„W (a) / dt W"^ (i) f (t) - W (x) D-^BfoW (6) / dt W"^ (i) f (i) , 

J a J a; 



or, 



where Green's matrix is given as 



Gi {x,t) = { 



\{x)= j dt Gi {x,t){{t) 

J a 

-W(x)D-iB6W(6) (t), x<t, 
W(a;)D-iBaW(a)W-i(i), x > t. 



(21) 



G{x,t) = < 



(22) 



2. Green's matrix of the second order system 
Using (16), Green's matrix for the second order system is 

-U„x2nW (x) D-lBfcW (6) (t) L2nxnP-' (t) , X< t, 

U„x2nW (x) D-iB,W (a) W-i (t) P-I(i), X>i. 

Note that W which is the fundamental matrix of the first order set is the Wronskian matrix of the 
second order set. With this result, we actually achieve our goal which is to construct the Green's 
matrix of the second order, self-adjoint matrix differential operator from the linearly independent 
solutions of the corresponding differential equation set. However, we continue to work on this result 
in order to find a more compact form and to find forms in which properties of the Green's matrix 
are more transparent. 

First of all, let's transform Green's matrix to a form in which boundary conditions is explicit. In 
order to achieve this goal, let's assume that the 2n linearly independent solutions forming Wronskian 
matrix, W is chosen in such a way that n of them satisfy one boundary condition, and the other n 
satisfy the other boundary condition; i.e. 

M^u" (x) = o, u" (a) = o, 
M^v^ (x) = o, v'^ (6) = o, 



where a, /3 = 1, . . . , n. Then, let's choose W in the form 



W 



U V 

u' v 
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where U and V are n x n matrices whose columns are u" and v^, respectively. Using this form of 
W greatly simplifies the matrix multiplications given in (22) and yields the results 



D-^BftW (6) 



I O 

o o 



D-iB„W(a) 



o o 

O I 



Block inverse of Wronskian matrix can be given as (See Sect.V); 

u-i + u-i V (V - u'u-i V) u'u-i -u-i V (V - U'U-i V) 
- (V - U'U-1 V) ~^ u'u-1 (V - U'U-1 V) 



w 



where each matrix is a function of t. Note that U (t) is singular at x = t = a. Thus, x < t part of 
the Green's matrix should be written in such a way that x = a boundary value is given separately. 
Putting these results in (22) yields a form where boundary conditions are explicitly satisfied: 
( 



O, a = X < t, 

G{x,t) = l\3 (a;) {t) V {t) (V (t) - U' (t) U-i (t) V {t)) P-i (t) , a < X < t, 



(23) 



y{x) (V (t) - U' (i) u-i (t) V (t))-' p-i (t) , 

Or, after rearranging, one has a more symmetric form; 



x>t. 



G(x,t) = < 



O 



u (x) u-i (t) [p (t) (v it) v-1 (i) - u' (t) u-1 (t))] ^ 
V (x) v-1 (t) [p (t) (v (t) (t) - u' (t) u-1 (t))] 



o 



a = X <t, 
a < X <t, 
b> x>t, 
b = x>t. 



(24) 



With the help of (12), a final compact form of the Green's matrix is obtained as 

U(x) (C^)"W^(i), x<t, 



G{x,t) = < 
Or, writing in terms of the elements; 

Gij{x,t) = {C-')^^{ 



(25) 



V (x) C-^U^ (i) , 



x>t. 



uf{x)v'^{t), x<t, 

(x) < (t), X > t. 



(26) 



where there is summation on the Greek indices. 

In all of the above forms of the Green's matrix, some properties are explicit, while the others 
not. Now, let's investigate these properties. 
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3. Verifying the properties of the Green's matrix 

Let's show that the Green's matrix that we have constructed satisfies the properties listed in 
Section II: 

• It obviously satisfies the homogeneous boundary conditions in the forms starting with (23). 

• It's columns are formed by the solutions of the homogeneous differential equation. 

• It is continuous at x = which is explicit in (23) and (24). 

• In forms (25) and (26), symmetry property is explicit. Let's write G^(t,a;); 



{t,x) = { 



which is simply equal to G (x, t). 



V (x) C-^U^ (t) , t< X, 
U(x) (C-i)^V^(i), t>x, 



• Derivative property of the Green's matrix can be given in the matrix form as; 

t+e 



lim -^G (x, t) 
e-»o dx 



Using the Green's matrix form given in (24); 



= it) ■ (27) 

t-e 



V (t) v-1 {t) [p {t) (V (t) v-1 (t) - u' (t) u-1 (t))] 

- u (x) u-^ (t) [p it) (V {t) it) - u' it) it))] = p-1 it) . 

C. Construction with a Guess 

Now, let's try to construct the Green's matrix in the similar way as in the case of a single 
differential operator. A somewhat similar derivation was given in [9] for a specific case. Since 
the columns of the Green's matrix are solutions of the homogeneous differential equation. Green's 
matrix should have the following form in order to satisfy the boundary conditions; 



G{x,t) 



U (x) S{t), x<t, 
Y{x)T{t), x>t, 



where S (t) and T [t) are unknown matrices. In this form, each column is a linear combination 
of the linearly independent solutions. Surely, this form is the most general guess satisfying the 
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boundary conditions and the homogeneous differential equation set. After putting the proposal in 
the derivative property in (27), one obtains 

V (t) T (t) - U' (t) S {t) = P-^ {t) . 

Using (12), one has 

p-i = (v'v-^ - u'U"^) vc-^u^. 

Putting it in the derivative property (27) and using the symmetry of VC~''^U'^ yield, 

V'T - U'S = V'C-^U^ - U' (C-^)^ V^. 
Then, the unknown S and T matrices are 

s = (c-^)^v^, 

T = C~^U^. 
This result yields the same form given in (25). 

IV. DISCUSSION AND CONCLUSION 

In this paper, we constructed the Green's matrix of a second order, self-adjoint matrix differ- 
ential operator. This construction is useful especially in the numerical studies, since obtaining the 
linearly independent solutions of the corresponding homogeneous differential equation set is easy 
with linearly independent initial conditions. Just after obtaining the set of linearly independent so- 
lutions, one may directly use (22) and obtain the Green's matrix without considering the boundary 
behavior of the solutions. However, once the linearly independent solutions are redefined such that 
half of them satisfy one homogeneous boundary condition and the other half the other boundary, 
or are obtained directly in this way with proper boundary conditions, then the compact form of 
the Green's matrix in (25) can be used by calculating the constant matrix in (7). By taking this 
route, one may avoid numerical errors coming from correspondingly greater amount of matrix mul- 
tiplications and inversions. Also, in cases where solving the corresponding homogeneous differential 
equation set analytically is easy, construction of the Green's matrix by use of (25) may become as 
easy as other techniques. 

Extracting Green's matrix of a higher order differential equation set from the corresponding first 
order differential equation set is a useful technique which is not well known in physics literature. 
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Although we construct the Green's matrix of a second order, self-adjoint matrix differential operator 
by using the Green's matrix of the corresponding first order differential equation set due to its 
physical relevance, Green's matrix for any higher order linear (matrix) operator, either having self- 
adjointness property or not, can be extracted from the Green's matrix of the corresponding first 
order differential equation set. 

A final comment on the boundary conditions is that a differential equation set satisfying bound- 
ary conditions other than the homogeneous ones can be handled as in the case of the single dif- 
ferential operator. Construction of a Green's matrix satisfying boundary conditions other than 
homogeneous ones can be handled by the first method given in this paper. 



APPENDIX A: INVERSE OF A BLOCK MATRIX 



We reproduce the derivation given by Thornburg [19]. Let's have a 2n x 2n matrix in the form 

A B 
C D 

where A, B, C, D are n x n matrices. Inverse of this matrix can be determined by obtaining the 
block decomposition of this matrix. 

Inverse of the following matrix forms can be found easily: 



A O 


-1 


A-i O 




O B 


-1 


o c-^ 


O D 




O D-i 




C O 




B-1 O 



I B 


-1 


I -B 




I O 


-1 


I o 


O I 




O I 


1 


C I 




-C I 



If we block decompose a general matrix in such a way that the above forms appear, then taking 
inverse can be handled by using the above relations. In order to obtain the block decomposition, 
one can use the following equation 

AE + BF = G, 
CE + DF = H. 

This equation set can be solved for F by multiplying the first equation with — CA~^ and summing 
with the second one. These operations are equal to multiplying coefficient matrix with 

I O 

CA-i I 



A 


B 




E 




G 


C 


D 




F 




H 
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from left. Then, coefHcient matrix becomes 

A B 

O D-CA-^B 

where = D — CA^^B is called Schur complement of A. In order to solve equations in E, the 
second equation should be multiplied with — BS^^ and summed with the first equation. These 
operations are equal to multiplying the modified coefficient matrix with 

I -BS-/ 
O I 

from left. Afterwards, the coefficient matrix becomes 

A O 
O Sa 

Thus, 



I 


-BS-/ ' 




I 


O 




A 


B 




A 


O 


o 


I 




-CA-i 


I 




C 


D 




O 





which yields the block decomposed form 



A 


B 




I 


O 




I BS-/ 




A 


O 


C 


D 




CA-i 


I 




O I 




O 


Sa_ 



Then, the inverse of the 2n x 2n matrix can be found as 



A 


B 


-1 


A-i + A-iBS;^^CA-i -A-^BS^^ 


C 


D 




-S-/CA-^ S-/ 



Using the above result, the inverse of the Wronskian matrix is; 



U V 


-1 


u-1 + u-i V (V - 


U'U 


-iv)"^u'u-i 


-u-iy (v - u'u-^v) ^ 


U' v 




- (v - u'u- 






(v'-u'u-^v)"^ 



Note that linear independence of the solutions of the corresponding homogeneous differential 
equation set implies invertibility of the Wronskian matrix and its four elements. 
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